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A density dependent M3Y effective nucleon-nucleon (NN) interaction which was based on the 
G-matrix elements of the Reid- Elliott NN potential has been used to determine the incompressibity 
of infinite nuclear matter. The nuclear interaction potential obtained by folding in the density 
distribution functions of two interacting nuclei with this density dependent M3Y effective interaction 
had been shown earlier to provide excellent descriptions for medium and high energy a and heavy ion 
elastic scatterings as well as a and heavy cluster radioactivities. The density dependent parameters 
have been chosen to reproduce the saturation energy per nucleon and the saturation density of 
spin and isospin symmetric cold infinite nuclear matter. The result of such calculations for nuclear 
incompressibility using the density dependent M3Y effective interaction based on the G-matrix 
elements of Reid-Elliott NN potential predicts a value of about 300 MeV for nuclear incompressibility. 
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The compression modulus or the incompressibility of infinite nuclear matter K^, denned through the curvature 
of the energy per nucleon E/A=e of the cold spin and isospin symmetric infinite nuclear matter at the saturation 
density is of fundamental importance. Different types of nuclear equation of state (EOS) considered in the study 
of central heavy ion collisions [1] are normally distinguished by different values. Experimentally, a widely used 
I , method is the determination of the incompressibility from the observed giant monopole resonances (GMR) [2-4] . But 
■ experimentally extracted incompressibility values [4] from the GMR data still show the span of values no better 
I/"") \ than 1.7 x (200-300) MeV, and more precise data are needed to pin down the range. Evidence from the various 
C***~ ■ nuclear data and neutron star masses favor a high compression modulus Koo (approximately equals) 300 MeV [5]. An 
accurate determination of Koq is very important for the study of properties of nuclei (radii, masses, giant resonances 
etc.), supernova collapses, neutron stars, and heavy ion collisions. 

A well-defined effecttive nucleon-nucleon (NN) interaction in the nuclear medium is important not only for different 
structure models but also for the microscopic calculation of the nucleon-nucleus and nucleus-nucleus potentials used 
in the analysis of the nucleon- and heavy-ion scattering. Effective NN interaction can be best constructed from a 
sophisticated G-matrix calculation. This interaction has been derived by fitting its matrix elements in an oscillator 
basis to those elements of the G-matrix obtained with the Reid-Elliott soft-core NN interaction [6]. The ranges of 
the M3Y forces were chosen to ensure a long-range tail of the one-pion exchange potential as well as a short range 
repulsive part simulating the exchange of heavier mesons. Such an effective NN interaction has been shown to provide 
a more realistic shape of the scattering potentials of the nucleon or heavy ion optical potentials obtained by folding 
, in the density distribution functions of two interacting nuclei with the effective NN interaction [7] . 

The real part of the nuclear interaction potential obtained by folding in the density distribution functions of 
two interacting nuclei with the M3Y effective interaction with one nucleon exchange evaluated from first principles 
was shown to provide good descriptions for medium and high energy a and heavy ion elastic scatterings [8,9]. It was 
observed that the high energy a and heavy ion elastic scattering data were better fitted, requiring less renormalisation 
and better chi-square per degrees of freedom (x 2 /F), when folding model potentials obtained from M3Y effective 
interaction based on the G-matrix elements of Reid-Elliott NN potential were used instead of those obtained from 
M3Y effective interaction based on the G-matrix elements of Paris NN potential. The semirealistic explicit density 
dependence on the M3Y effective interaction with a zero-range pseudo-potential was then employed. The zero-range 
pseudo-potential represents the single-nucleon exchange term while the density dependence takes care off the higher 
order exchange effects and the Pauli blocking effects. Since the density dependence of the effective projectile-nucleon 
interaction has been found to be fairly independent of the projectile, as long as the projectile-nucleus interaction 
is amenable to a single-folding prescription, implies that in a double folding model, the density dependent effects 
on the nucleon-nucleon interaction can be factorized into a target term times a projectile term [10,11]. The folding 
model potentials obtained using such factorized density dependent M3Y-Reid-Elliott effective interaction along with 
a zero-range pseudo-potential were then used to describe successfully the medium and high energy a and heavy ion 
clastic scatterings [12,13]. 
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Subsequently, the folding model potentials obtained using M3Y effective interaction based on the G-matrix elements 
of Reid-Elliott NN potential were then used to describe successfully the heavy particle decay [14] and a decay [15] 
while it was found to provide excellent description for cluster radioactivity when the spin-parity conservation was 
also taken into account [16]. In the present work nuclear incomprcssibility has been calculated theoretically using 
the same density dependent M3Y effective interaction, based on the G-matrix elements of Reid-Elliott NN potential, 
supplemented by a zero range pseudo-potential. It is worthwhile to mention here that due to attractive character of 
the M3Y forces the saturation condition for cold nuclear matter is not fulfilled. However, the realistic description 
of nuclear matter properties can be obtained with the density dependent M3Y effective interaction. Therefore, the 
density dependent parameters have been obtained by reproducing the saturation energy per nucleon and the saturation 
nucleonic density of the spin and isospin symmetric cold infinite nuclear matter. 

The incomprcssibility of the spin and isospin symmetric cold infinite nuclear matter is defined as 

K x = k 2 F d 2 e/0k 2 F = 9p 2 d 2 e/dp 2 \ p=Pao (1) 

where the Fermi momentum k F for the spin and isospin symmetric infinite nuclear matter is given by 

k% = 1.5tt 2 p (2) 

and p is the nucleonic density while p^ being the saturation density for the spin and isospin symmetric cold infinite 
nuclear matter. 

The effective nucleon-nucleon interaction v(s) is assumed to be density and energy dependent and therefore becomes 
functions of density and energy and is generally written as [12,13] 

v(s,p,e)=t M3Y ( S ,e)g(p,e) (3) 

where t M3Y is given by [6] 

t M3Y (s,e) = 7999exp(-4s)/(4s) - 2134exp(-2.5s)/(2.5s) + J m (e)5(s) (4) 

where the zero-range pseudo-potential representing the single- nucleon exchange term is given by [17] 

Joo(e) = -276(1 - 0.005e)(MeV./m 3 ) (5) 

and the density dependent part has been taken to be [18] 

g(p,e) = C(l-(3(e)p 2 / 3 ) (6) 

which takes care of the higher order exchange effects and the Pauli blocking effects. 

The energy per nucleon e obtained using the effective nucleon-nucleon interaction v(s) for the spin and isospin 
symmetric cold infinite nuclear matter is given by 

e = 3h 2 k 2 F /10m + g(p,e)pJ v /2 (7) 

where m is the nucleonic mass equal to 931.4943 MeV/c 2 and J v represents the volume integral of the M3Y interaction 
supplemented by the zero-range pseudopotential having the form 

J v (e) = J J y"t M3y (s,e)d 3 s = 7999(47r/4 3 )- 2134(47r/2.5 3 ) + Joo(e) (8) 

The equilibrium density of the nuclear matter is determined from the saturation condition de/dp = which yields 
the equation 

de/dp = h 2 k 2 F /5mp + J V C(1 - 5(3(E)p 2 / 3 /3)/2 = (9) 

while the eqn.(7) can be rewritten with the help of eqn.(6) as 

3h 2 k 2 F /Wm + pJ v C(l - (3(E)p 2 ' 3 )/2 = e (10) 

so that eqn.(9) and eqn.(10) can be solved simultaneously for any values of the saturation energy per nucleon and the 
saturation density of the spin and isospin symmetric cold infinite nuclear matter to obtain the values of the density 
dependent parameters (3(e) and C. Density dependent parameters (3(e) and C, thus obtained, can be given by 
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/?(e) = [(3-3p)/(9-5p)]/p 2 /3 



where 

p= [10me]/[?i 2 (1.57r 2 p) 2/3 ] (12) 

and 

C = -[2h 2 k 2 F ]/[5mJ vP (l - 5/3(e)/9 2 / 3 /3)] (13) 

respectively. 

Once the density dependent parameters /3(e) and C have been fixed for a saturation energy per nucleon and a 
saturation density of a spin and isospin symmetric cold infinite nuclear matter, the incompressibility can be 
evaluated from the equation given by 

K x = [-3h 2 k 2 F /5m - bJ v CI3(e) P b '\ =Poo (14) 

which has been obtained using eqn.(l), eqn.(2) and eqn.(9). 

The results of the present calculations using the saturation energy per nucleon e = -16 MeV for the spin and 
isospin symmetric cold infinite nuclear matter have been presented in Table- 1. Since the value used for the saturation 
density by different groups do differ, a narrow range of its acceptable values have been used. The value of the density 
dependence parameter (3(e), which is supposed to dependend on energy, has been found rather close to 1.6 fm 2 , a 
value which was obtained by optimum fit to the experimental data for a radioactivity which involves very low energies 
[16] compared to the high energy a elastic scattering [12]. Using the usual value of 0.005/MeV for the parameter 
of energy dependence of the zero range pseudo-potential, as one can see in eqn.(5), the value obtained for the other 
density dependence parameter C is close to 2. However, readjusting the value of the energy dependence parameter to 
0.100/MeV brings down the value of the density dependence parameter C approximately equal to 1, while the value of 
the nuclear incomressibility remains unaltered which is obvious from eqn.(13) and eqn.(14). The values obtained for 
the nuclear incompressibility ranges from 304 MeV - 310 MeV for acceptable range of values of saturation densities 
used here. 



TABLE I. Compression modulus at different saturation densities 



Poo 


fm 2 


C 


Energy dependence parameter 




fm- 3 




MeV' 1 


MeV 


0.170 


1.551 


1.98 


0.005 


309.6 


0.170 


1.551 


1.02 


0.100 


309.6 


0.165 


1.586 


2.02 


0.005 


308.2 


0.165 


1.586 


1.04 


0.100 


308.2 


0.160 


1.624 


2.07 


0.005 


306.9 


0.160 


1.624 


1.06 


0.100 


306.9 


0.155 


1.664 


2.11 


0.005 


305.5 


0.155 


1.664 


1.09 


0.100 


305.5 


0.150 


1.705 


2.16 


0.005 


304.0 


0.150 


1.705 


1.11 


0.100 


304.0 
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Theoretical estimate of Koo from refractive a-nucleus scattering claims to have pinned down its range to about 240 
MeV-270 MeV [19,20] ruling out the lower value of about 180 MeV. Another recent theoretical estimate by infinite 
nuclear matter model (INM) [21] claims a well defined and stable value of = 288 ± 20 MeV. The finite nuclear 
incompressibility Ka can be related to the energies Eqmr of the giant monopole resonance (GMR) as 

Ka — M < r 2 > E 2 GMR /h 2 (15) 

where M is the mass of the nucleus, < r 2 > is its mean square radius. Following a liquid drop model type expansion, 
one writes down a similar expression for Ka as 

K A = K V + KsA- 1 / 3 + K CV A- 2 ' Z + K asym (l - 2Z/A) 2 + K C Z 2 A^ (16) 

and then the volume term K v can be identified as K^. Although eqn.(16) looks to be a natural path to obtain 
from Ka, this expansion is beset with many difficulties in regard to its convergence and the number of data on GMR 
being limited, no precise value of if oo can be extracted [4] . A recent determination of Koo based upon the production 
of hard photons on heavy ion collision led to the experimental estimate K^ = 290 ± 50 MeV [22] and the result of 
the present calculation for K^ is in excellent agreement with this recent experimental finding. 

A density dependent M3Y-Reid-Elliott effective nucleon-nucleon interaction which has been used within a folding 
model prescription to provide excellent description for medium and high energy a and heavy ion scatterings, a and 
heavy ion cluster radioactivities, has been used to calculate the nuclear incompressibility K^ of the spin and isospin 
symmetric cold infinite nuclear matter. The resulting value of K^ is found to be 304 MeV-310 MeV approximately. 
This value is in excellent agreement with the recent experimental estimate of K^ = 290 ± 50 MeV and rules out lower 
values of about 180 MeV. Results of the present calculations are very close to the recent theoretical estimate by INM 
of 288 ± 20 MeV. 
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